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Abstract 

We present the applications of variational-wavelet ap- 
proach for computing multiresolution/multiscale represen- 
tation for solution of some approximations of Vlasov- 
Maxwell equations. 



level/scale of resolution. Starting in part 2 from Vlasov- 
Maxwell equations we consider in part 3 the generaliza- 
tion of our approach based on variational formulation in 
the biorthogonal bases of compactly supported wavelets. 



1 INTRODUCTION 

In this paper we consider the applications of a new numeri- 
cal-analytical technique which is based on the methods of 
local nonlinear harmonic analysis or wavelet analysis to the 
nonlinear beam/accelerator physics problems described by 
some forms of Vlasov-Maxwell (Poisson) equations. Such 
approach may be useful in all models in which it is possi- 
ble and reasonable to reduce all complicated problems re- 
lated with statistical distributions to the problems described 
by systems of nonlinear ordinary/partial differential equa- 
tions. Wavelet analysis is a relatively novel set of math- 
ematical methods, which gives us the possibility to work 
with well-localized bases in functional spaces and gives for 
the general type of operators (differential, integral, pseu- 
dodifferential) in such bases the maximum sparse forms. 
Our approach in this paper is based on the generalization 
of variational-wavelet approach from [l]-[8], which allows 
us to consider not only polynomial but rational type of non- 
linearities [9]. The solution has the following form (related 
forms in part 3) 
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which corresponds to the full multiresolution expansion in 
all time/space scales. 

Formula ([]]) gives us expansion into the slow part u 
and fast oscillating parts for arbitrary N. So, we may move 
from coarse scales of resolution to the finest one for ob- 
taining more detailed information about our dynamical pro- 
cess. The first term in the RHS of formulae (1) corresponds 
on the global level of function space decomposition to reso- 
lution space and the second one to detail space. In this way 
we give contribution to our full solution from each scale 
of resolution or each time/space scale. The same is correct 
for the contribution to power spectral density (energy spec- 
trum): we can take into account contributions from each 
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Figure 1: Base wavelet. 



2 VLASOV-MAXWELL EQUATIONS 

Analysis based on the non-linear Vlasov-Maxwell equa- 
tions leds to more clear understanding of the collecti- 
ve effects and nonlinear beam dynamics of high inten- 
sity beam propagation in periodic-focusing and uniform- 
focusing transport systems. We consider the following 
form of equations ([11] for setup and designation): 
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The corresponding Hamiltonian for transverse single-par- 
ticle motion is given by 



1 



H(x,y,x',y',s) = -(x' 2 +y> 2 ) 
+\[kx{s)x 2 + k y (s)y 2 } + if)(x, y, s). 



(4) 



Related Vlasov system describes longitudinal dynamics 
of high energy stored beam [12]: 
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I(6;T)= f dvvf(0,v;T). 
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3 VARIATIONAL APPROACH IN 
BIORTHOGONAL WAVELET BASES 

Now we consider some useful generalization of our varia- 
tional wavelet approach. Because integrand of variational 
functionals is represented by bilinear form (scalar product) 
it seems more reasonable to consider wavelet constructions 
[13] which take into account all advantages of this struc- 
ture. The action functional for loops in the phase space is 



F( 7 ) 



pdq 



H(t,j(t))dt 



(8) 



The critical points of F are those loops 7, which solve the 
Hamiltonian equations associated with the Hamiltonian H 
and hence are periodic orbits. Let us consider the loop 
space ft = C oc (S 1 ,R 2n ), where S 1 = R/Z, of smooth 
loops in R 2n . Let us define a function $ : — > R by 
setting 

= ( 7, < - J x> x > dt- { H(x(t))dt, x eil 
Jo 2 Jo 

(9) 

Computing the derivative at x <G ft in the direction of y e 
il, we find 
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< —Jx - SjH(x), y > dt 



(10) 



Consequently, $'(x)(y) = for all y e iff the loop x is 
a solution of the Hamiltonian equations. Now we need to 
take into account underlying bilinear structure via wavelets. 
We started with two hierarchical sequences of approxima- 
tions spaces [13]: ■■■V-2 C V-i C V C V\ C V 2 ■ ■ ■ , 
. . . V-2 C V-i C V C V\ C V 2 . . . , and as usu- 
ally, W is complement to Vo m Vi> t> ut now not neces- 
sarily orthogonal complement. New orthogonality condi- 
tions have now the following form: Wq _L Vq, Wq _L 



_L Wj , V, _L W 3 translates of tp span Wq, 



translates of tp span Wo- Biorthogonality conditions 
are < tpj k ,^j'k' >= $- 00 i>jk{x)'ipj'k'{x)&x = <5 fc fe'%', 
where tpjk (x) = 2i/ 2 tp(2 : > x— k). Functions <p(ar), (p{x — k) 
form dual pair: < <^(a; — k),ip(x — I) >= Ski, < 
(f(x — k),ip(x — £) >= for Vfc, W. Functions ip,ip 
generate a multiresolution analysis. ip(x — k), ip(x — k) 
are synthesis functions, <p(x — i), tp(x — I) are analysis 
functions. Synthesis functions are biorthogonal to analysis 
functions. Scaling spaces are orthogonal to dual wavelet 
spaces. Two multiresolutions are intertwining Vj + Wj = 
Vj+\, Vj + Wj = Vj+i. These are direct sums but not 
orthogonal sums. So, our representation for solution has 
now the form 



/(*) = EiwM*), 



where synthesis wavelets are used to synthesize the func- 
tion. But bjk come from inner products with analysis 



wavelets. Biorthogonality yields 

b lm = f f(t)$ tm (t)dt. (12) 

So, now we can introduce this more useful construction 
into our variational approach. We have modification only 
on the level of computing coefficients of reduced nonlinear 
algebraical system. This new construction is more flexi- 
ble. Biorthogonal point of view is more stable under the ac- 
tion of large class of operators while orthogonal (one scale 
for multiresolution) is fragile, all computations are much 
more simpler and we accelerate the rate of convergence. In 
all types of (Hamiltonian) calculation, which are based on 
some bilinear structures (symplectic or Poissonian struc- 
tures, bilinear form of integrand in variational integral) this 
framework leads to greater success. 

So, we try to use wavelet bases with their good spatial 
and scale-wavenumber localization properties to explore 
the dynamics of coherent structures in spatially-extended, 
'turbulent'/stochastic systems. After some ansatzes and re- 
ductions we arrive from (2), (3) or (5)-(7) to some system of 
nonlinear partial differential equations [10]. We consider 
application of our technique to Kuramoto-Sivashiinsky 
equation as a model with rich spatio-temporal behaviour 
[14] (0 < x < L, £ = x/L, u(0,t) = u(L,t), 
u x (0,t) = u x (L,t)): 
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where a\{€) are both wavelets. 

Variational formulation 
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reduces (13) to ODE and algebraical one. 
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(11) In particular case on V2 \ Vq we have: 





Figure 2: The solution of eq.(13) 



Figure 3: The solution of eq.(17) 



ca^CLi — caoa2 + da\ — da 2 . 
— coq — daoai + £aoa,2 — fai<i2 ~~ f a 2 
ca 2 — lagai + da^ai + da^a2 + ja\ + fa\a^ 

Then in contrast to [14] we apply to (16) methods from [1]- 
[9] and arrive to formula (1). The same approach we use 
for the general nonlinear wave equation 



where 



u tt = Uxx ~ mu - f(u), 



f(u) =au 3 + J2 hu k 

k>5 



(17) 



(18) 



According to [2], [10] we may consider it as infinite dimen- 
sional Hamiltonian systems with phase space = Hq x L 2 
on [0, L] and coordinates: u, v = u t , then 

ii r 

H = — < v,v > +— < Au,u > + / g(u)dx 
2 2 Jo ' 

d 2 f 
A = + m, g= I /(s)ds 

u t = -5-=« (19) 

ov 

dH a *r \ 
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or u(t) = JVK(u(t)) 
Then anzatzes: 

u(t,x) = U(ujit, . . . ,uj n t,x) (20) 

u(t,x) = Uk{x) exp(ifc • uj(k)t) 

kez n 

u(t, x) = S(x — vt) 
u(t,x) = Uk{x)V k {t) 

and methods [ 1 ] - [ 1 0] led to formulae (1). Resulting mul- 
tiresolution/multiscale representation in the high-localized 
bases (Fig.l) is demonstrated on Fig. 2, Fig. 3. We would 
like to thank Professor James B. Rosenzweig and Mrs. 
Melinda Laraneta for nice hospitality, help and support dur- 
ing UCLA ICFA Workshop. 
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